Structural and magnetic phase transitions in Nai.^FeAs 
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We use neutron scattering to study the spin and lattice structures of single crystal and powder 
samples of Nai-jFeAs (Tc = 23 K). On cooling from room temperature, the system goes through 
a series of phase transitions: first changing the crystal symmetry from tetragonal to orthorhombic 
at 49 K, then ordering antiferromagnetically with a spin structure similar to that of LaFeAsO and 
a small moment (0.09±0.04 /is), and finally becoming superconducting below about 23 K. These 
results confirm that antiferromagnetic order is ubiquitous for the parent compounds of the iron 
arsenide superconductors, and suggest that the separated structural and magnetic phase transition 
temperatures are due to the reduction in the c-axis exchange coupling of the system. 



Determining the universal features of iron arsenide su- 
perconductors is an important first step in developing 
a microscopic theory to understand the high-transition 
temperature (high- Tc) superconductivity in these mate- 
rials From the outset, it was known that antifer- 
romagnetic (AF) order is a competing ground state to 
superconductivity in iron arsenide superconductors. The 
parent compound LaFeAsO exhibits a structural distor- 
tion at 155 K and then orders antiferromagnetically be- 
low 137 K, electron doping to induce superconductivity 
suppresses both the structural distortion and static AF 
order 2, 3] • Although subsequent neutron scattering ex- 
periments on the parent compounds of other iron-based 
superconductors including i?FeAsO (i? = Ce, Pr, Nd) 
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^FesAsa {A = Ba, Sr, Ca) [1, i, [3, and 
have found similar lattice distortion and AF 
order as that of the LaFeAsO, the MFeAs system (M — 
Li, Na) seemed to be an exception to this universal pic- 
ture since the initial neutron and X-ray scattering ex- 
periments have found no evidence of lattice distortion or 
static AF order [H [M [H, [11]. These results are in 
contrast with local density approximation (LDA) calcu- 
lations, where the Fermi surfaces and magnetic orders of 
MFeAs are expected to be similar to that of the LaFeAsO 
[TtI . [ist . Although recent transport and heat capacity 
measurements on single crystals of Nai-^FeAs showed 
two anomalies at 52 K and 41 K that are assigned to 
structural and AF phase transitions, respectively [l^, 
muon-spin rotation experiments (/LtSR) confirmed only 
the AF ordering and neutron scattering found no evi- 
dence for structural distortion 

In this paper, we report neutron scattering investiga- 
tion of spin and lattice structures of single crystals and 
polycrystalline materials of Nai-^FeAs. We identify a 
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tetragonal to orthorhombic structural phase transition 
near 50 K and thus confirm the transport measurements 
[l^ . Although our neutron powder diffraction measure- 
ments were unable to detect AF order due to the small Fe 
moment, single crystal experiments using thermal triple- 
axis spectroscopy unambiguously confirmed an AF phase 
transition below 40 K and showed that the system forms 
a coUinear in-plane AF spin structure identical to other 
iron arsenides [1, [1, [1, [3,11, i, ^ but doubling the unit 
cell along the c-axis [Fig. 1(a)]. The ordered moment is 
by far the smallest in iron arsenides, being 0.09 ± 0.04 /ib- 
These results suggest that AF order with a collinear in- 
plane spin structure is ubiquitous property of the parent 
compounds of iron arsenide superconductors. We argue 
that the separated structural and magnetic phase tran- 
sitions in Nai-^FeAs is due to reduced c-axis exchange 
coupling and Na deficiencies. Bulk superconductivity in 
Lii-^FeAs and Nai-^FeAs can only arise with enough 
self-doping induced by alkali metal deficiencies. 

We prepared 7 grams of polycrystalline Nai-^FeAs 
sample and 0.6 gram of single crystals as described in 
Ref. [ll]. The resistivity measurement gives the on- 
set and zero- resistivity Tc as 23 K and 8 K, respec- 
tively ^1^. Powder neutron diffraction measurements 
were performed on the BT-1 high resolution powder 
diffractometer at the NIST Center for Neutron Research. 
The BT-1 diffractometer has a Ge(3,l,l) monochroma- 
tor (A = 2.0785 A) and collimators with horizontal 
divergences of 15'-20'-7' full width at half maximum 
(FWHM). Powder diffraction data refinement was done 
by using the GSAS program. The measurements on 
single crystals were carried out on the HB-l triple- axis 
spectrometer at the High Flux Isotope Reactor, Oak 
Ridge National Laboratory. We fixed final neutron en- 
ergy at Ef — 13.5 meV and used PG(0,0,2) (pyrolytic 
graphite) as monochromator and analyzer. The coUi- 
mations for magnetic and structural measurements are 
48'-60'-80'-120' and 15'-20'-20'-30', respectively. We de- 
note Q(A-i) = {2nH/a,2nK/b,2TTL/c), where ar = 
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TABLE I: Refinement results of powder diffraction data 
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FIG. 1: (a) Nuclear and (b) magnetic structures of the ideal 
NaFeAs. (a) includes two orthorhombic nuclear unit cells for 
comparison with the magnetic unit cell in (b). (c) and (d) 
show the [H, 0, 1] and [1,0,1/] scans around (1,0, 1)m mag- 
netic Bragg peak at 5 K and 70 K. Both peaks disappear at 
70 K, indicating their magnetic nature. 



br = 3.94481(3) A, ct = 6.99680(8) A for tetragonal 
structure and ao = 5.58906(8) A, bo = 5.56946(8) A, 
Co = 6.9919(1) A for orthorhombic structure. The mag- 
netic unit cell is defined as ao x bo x ^co- 

We first describe neutron powder diffraction measure- 
ments. At a temperature well above the structural and 
magnetic phase transitions (T = 70 K) , Rietvelt analysis 
reveals a tetragonal structure with space group PA/nmm 
consistent with earlier results After cooling the 

sample down to 5 K, the tetragonal (2, 2, 1) peak splits 
into two peaks as shown in Fig. 2(a). Detailed tem- 
perature dependence of the (2,2, 1) profiles in Fig. 2(a) 
reveals that the tetragonal to orthorhombic structural 
phase transition occurs near 45 K. Refinement of the 
diffraction pattern supports an orthorhombic structure 
at low temperature and detailed structural parameters 
are listed in Table 1 for the two temperatures investi- 
gated. Fixing the occupancies of Fe and As to 1, we 
obtain the Na content as 0.985(7), or equivalently 1.5% 
Na deficiencies. Figures 2(b)-2(c) show the temperature 
dependence of some key parameters. As a function of 
increasing temperature, the Fe-Fe bond length decreases 
until they become equal (tetragonal) while the Fe-As dis- 
tance remains unchanged [Fig. 1(c)]. The orthorhombic 
lattice parameters ao and bo behave similarly [Fig. 1(b)]. 
While the nearest neighbor bond angles change in oppo- 
site directions with increasing temperature, the diagonal 
bond angle is essentially temperature independent [Fig. 
2(d)]. These results are similar to those of LaFeAsO 21|, 
and thus suggesting the same underlying mechanism for 



Nao.985FeAs(5 K), Cmma, x = 1-453 
ao = 5.58906(8) A, bo = 5.56946(8) A, cq 



6.9919(1) A 



Atom site 

Na 4a 

Fe 4g 

As 4a 



X 


0.25 




y z 

0.25 0.3533(2) 



0.25 0.7977(1) 



occupancy 
0.985(7)" 
1^ 



Selected bond lengths and angles: 

Fe-Fe X 2 2.79453(4) A Fe-Fe x 2 2.78473(4) A 



Fe-Asx4 2.4272(4) A 
Na-Asx2 2.9874(6) A 
Fe-As-Fex2 70.29(1)° 
Fe-As-Fex2 108.72(3)° 



Na-Asxl 3.107(2) A 
Na-Asx2 2.9782(6) 
Fe-As-Fex2 70.01(1)° 



Nao.985FeAs(70 K), PA/nmm, x = 1-685 



At = bT = 


3.94481(3) A, 


Ct = 


6.99680(8) A 


Atom site 


X 


y 


z occupancy 


Na 2c 


0.25 


0.25 


0.3535(2) 


Fe 2a 


0.75 


0.25 


l' 


As 2c 


0.25 


0.25 


0.7976(1) l'' 



Selected bond lengths and angles: 
Fe-Fex4 2.78941(3) A Fe-Asx4 2.4282(4) A 
Na-Asx4 2.9830(6) A Na-Asxl 3.107(2) A 
Fe-As-Fex4 70.11(2)° Fe-As-Fex2 108.64(3)° 

"Na occupancy is calculated as the mean value of those at several 
temperatures. 
'Fe and As occupancies are fixed to 1. 
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FIG. 2: (a) Temperature dependence of the 29 scans around 
the (2, 2, 1)t Bragg peak in powder diffraction data across the 
tetragonal to orthorhombic phase transition. Clear splitting 
of the peaks is seen below 33 K. We collected full powder 
diffraction spectra at several temperatures and refinement re- 
sults gave the temperature dependence of several key param- 
eters: (b) lattice constants a, b and c; (c) bond lengths of 
Fe-Fe and Fe-As; (d) bond angles of Fe-As-Fe. In all panels, 
left and right axes have the same scale. 
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FIG. 3: (a) Temperature dependence of the FWHM of 61- 
20-scan at the nuclear peak (1,1, 0)r position using A/2 scat- 
tering by removing the pyrolytic graphite (PG) filter. The 
peak width clearly increases below about 50 K as marked by 
the vertical line, (b) Temperature dependence of the peak 
intensity at the AF peak (1, 0, 3) m position suggesting a Neel 
temperature of 39 K. The black line is the fitted result for 
the whole temperature range by using a simple Ising model 
described in the text, while the blue line only focuses on tem- 
peratures above 20 K. No anomaly is seen across Tc, sug- 
gesting that superconductivity is filamentary and not a bulk 
phenomenon. 



the structural phase transition. 

To precisely determine the structural transition tem- 
perature, we carefully measured the temperature depen- 
dence of the (1,1, 0)t peak width (in FWHM) on the 
single crystal sample using A/2 as shown in Fig. 3(a). 
Although the resolution of the triple-axis spectrometer 
is not good enough to resolve two separate peaks from 
the 9 — 29 scans at low temperature, the change of the 
FWHM reveals a structural phase transition temperature 
of To ~ 50 K, which is consistent with the higher tran- 
sition temperature in the transport measurements [loj. 

To search for possible magnetic order, we carried out 
measurements using triple-axis spectroscopy on both the 
powder samples and single crystals. While we cannot 
find any magnetic peak in the powder diffraction data 
due to small Fe moment, we observe clear AF order on 
single crystals at low temperatures. It turns out that 
the in-plane AF unit cell of Nai-^FeAs is identical to 
that of LaFeAsO 0, where the magnetic Bragg peaks 
are observed at (1,0, L)m {L — 1,3,5,7). Figures 1(c) 
and 1(d) show wave vector scans along the orthorhombic 
[H, 0, 1] and [1, 0, L] directions [6^ at 5 K and 70 K. The 



resolution limited peaks around (1,0, 1)m at 5 K disap- 
pear at 70 K [Figs. 1(c) and 1(d)]. Figure 3(b) shows 
the temperature dependence of the scattering at wave 
vector Q — (1,0, 3)mj where we estimate that the onset 
magnetic transition temperature is about 37 K with less 
than 1 K thermal hysteresis [Fig. 3(b)]. These results 
are consistent with transport measurements where the 
^--^40 K transition is identified as magnetic in nature. In 
a simple Ising model, the magnetic order parameter is 
related to temperature via ^(r)^ oc (1 - T/rAr)^^. For 
BaFe2As2, AF phase transition has a critical exponent 
/3 = 0.103 ± 0.018 0. Fitting the whole temperature 
dependence of (1,0,3)m peak intensity yields an unrea- 
sonable T/v = 34.7±0.9 K, shown as the black line in Fig. 
3(b). Limiting the fitting range to temperatures above 
20 K gives a Tat = 37.1 ± 0.2 K and /3 = 0.28 ± 0.02, 
considerably higher than that of BaFe2As2. However, 
we caution that the order parameter was measured using 
peak intensity on single crystals, rather than integrated 
intensity measurement. 

A collinear AF structure could either have spin direc- 
tions along the a- or 6-axis, which would correspond to 
magnetic peaks at (1, 0, L)m or (0, 1, L)m directions, re- 
spectively. In previous work on SrFe2As2 [ol, the AF 
ordering direction was determined to be along the a-axis 
direction by comparing the (1, 0, 1) magnetic bragg peak 
with A/2 scattering from (2,0,1) and (0,2,1) nuclear 
Bragg peaks. Since the orthorhombic peak splitting in 
the case of Nai-^FeAs is rather small, we used very tight 
collimations for this purpose. Figure 4(a) shows the mag- 
netic Bragg peak near (1,0, L)m or {0,1, L)m and A/2 
scattering from twinned peaks of (1, 0, 0.5) and (0, 1, 0.5). 
It is immediately clear that the magnetic (1,0, 1)m peak 
position in 6* — 26* scan is smaller than that of the nu- 
clear peak. This result is consistent with earlier work on 
SrFe2As2 i- 

To determine the spin direction, we calculate the mag- 
netic structure factors by assuming that the moments 
point to the a-axis direction. The observed magnetic in- 
tensities are obtained by integrating the 9-29 scans at the 
expected magnetic peak positions (1, 0, L) {L = 1, 3, 5, 7) 
in the three-axis mode. A comparison of the calculated 
and observed magnetic peak intensities shown in Fig. 
4(b) reveals that such a model explains the data reason- 
ably well. The small deviation between the observed and 
calculated intensities sets a limit of the moment direction 
to be within 15 degrees away from the a-axis. Therefore, 
the magnetic structure in Nai-^FeAs is the same as that 
in the AFe2As2 system as shown in Fig. 1(b) [sl. [ol. [l0| . 
Assuming this spin structure, we can estimate an Fe mo- 
ment of 0.09±0.04 fj,B by comparing intensities of the AF 
Bragg peaks with with a series of nuclear peaks. Such a 
small moment explains why we as well as another group 
[2(| cannot find any AF Bragg peaks in neutron powder 
diffraction measurements. 

Our results on Nai-^FeAs suggest that orthorhombic 
structure and collinear AF spin ordering in Fig. 1(b) 
are ubiquitous properties of undoped iron arsenides. Be- 
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FIG. 4; (a) Comparison of the 9-26 scans at (1,0, 1)m with 
and without PG fiher, measuring magnetic and nuclear peaks, 
respectively, (b) Calculated and experimental intensities of 
magnetic peaks (1, 0, L)m {L = 1, 3, 5, 7). In the calculation, 
the moment is assumed to be along ao axis. 



cause of the difficulty in making stoichiometric samples 
of NaFeAs, Na deficiencies in the as-grown NaFeAs dopes 
holes onto the FeAs-plane that can induce superconduc- 
tivity (IS] , For BaFe2As2, transport and neutron scatter- 
ing experiments have shown that electron-doping reduces 
the c-axis magnetic exchange coupling and separates the 
structural and magnetic phase transitions 23, 24, 25, 26]. 
Theoretically, it has been argued that the strength of 



the c-axis magnetic coupling controls the simultaneous 
or separated structural/magnetic phase transitions [27| . 
For lightly doped BaFei.96Nio.o4As2, inelastic neutron 
scattering experiments showed a dramatic drop in the 
c-axis correlations with electron-doping. Based on den- 
sity function theory calculations, the c-axis exchange 
coupling of NaFeAs is smaller than that of BaFe2As2 
but larger than that of LaFeAsO. Experimentally, the 
temperature separations between structural and mag- 
netic phase transitions are similar for LaFeAsO 0] and 
NaFeAs. In addition, the AF order parameter showed no 
anomaly across Tc similar to those of BaFei.96Nio.o4As2 
[Fig. 3(b)]. These results suggest that superconduc- 
tivity in our Nai-^FeAs is filamentary and not a bulk 
phenomenon. Since our single crystals of Nai-^FeAs is 
slightly doped away from ideal stoichiometry (Na defi- 
ciency), it is unclear whether the observed small moment 
and large differences in structural/magnetic phase tran- 
sition temperatures are the intrinsic or doping-induced 
c-axis coupling reduction. Future inelastic neutron scat- 
tering experiments will be able to determine the exchange 
coupling along the c-axis. 

In conclusion, we have determined the lattice and mag- 
netic structures of single crystal Nai-^FeAs. Our re- 
sults indicate that the parent materials of NaFeAs and 
LiFeAs superconductors have orthorhombic lattice dis- 
tortion and coUinear AF order. This work establishes 
that the orthorhombic structure and AF coUinear order 
are ubiquitous to all undoped iron arsenide materials. 
Superconductivity arises from electron or hole-doping of 
their AF parent compounds, and therefore suggesting 
that spin fluctuations are important for superconductiv- 
ity of these materials. 

We thank Jiangping Hu, J. A. Fernandez-Baca, Tao 
Xiang, Zhong-Yi Lu for helpful discussions. This work 
is supported by the U.S. NSF No. DMR-0756568, U.S. 
DOE BES No. DE-FG02-05ER46202, and by the U.S. 
DOE, Division of Scientific User Facilities. The work in 
lOP is supported by the Ministry of Science and Tech- 
nology of China and NSFC. 



[1] Y. Kamihara et al, J. Am. Chem. See. 130, 3296 (2008). 

[2] C. de la Cruz et al, Nature (London) 453, 899 (2008). 

[3] M. A. McGuire et al, Phys. Rev. B 78, 094517 (2008). 

[4] J. Zhao et al, Nature Materials 7, 953 (2008). 

[5] S. A. J. Kimber et al, Phys. Rev. B 78, 140503 (2008). 

[6] J. Zhao et al, Phys. Rev. B 78, 132504 (2008). 

[7] Y. Chen et al., Phys. Rev. B 78, 064515 (2008). 

[8] Q. Huang et al, Phys. Rev. Lett. 101, 257003 (2008). 

[9] J. Zhao et al, Phys. Rev. B 78, 140504(R) (2008). 
[10] A. 1. Goldman et al, Phys. Rev. B 78, 100506(R) (2008). 
[11] W. Bao et al, arXiv: 0809.2058. 
[12] ShiUang Li et al, Phys. Rev. B 79, 054503 (2009). 
[13] X. C. Wang et al, SoUd State Commun. 148, 538 (2008). 
[14] M. J. Pitcher et al. Chemical Communications, 5918 
(2008). 

[15] J. H. Tapp et al, Phys. Rev. B 78, 060505(R) (2008). 



[16] C. W. Chu et al, arXiv:0902.0806 
[17] D. J. Singh, Phys. Rev. B 78, 09451 1 (2008). 
[18] R. A. Jishi and H. M. Alyahyaei, ar Xiv:0812.1215l f2008'). 
[19] G. F. Chen et al, arXiv:0902.1100vl. 
[20] D. R. Parker et al. Chemical Communications, 2189 
(2009). 

[21] T. Nomura et al, Supercond. Sci. Technol. 21, 125028 

(2008). 

[22] S. D. Wilson et g^.. larXiv: 0903.2430Vl. 

[23] N. Ni et al, Phys. Rev. B 78, 214515 (2008). 

[24] J.-H. Chu et al, Phys. Rev. B 79, 014506 (2009). 

[25] C. Lester et al, Phys. Rev. B 79, 144523 (2009). 

[26] L. W. Harriger et al, arXiv:0904.3775vl. 

[27] C. Fang et al, Phys. Rev. B 77, 224509 (2008). 



